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(D ' Abstract 

Several years ago, Stefan Pokorski, Manfred Miinz and us outlined a program for calculation 

* I — — 

^ , of the NLO QCD corrections to the weak radiative B meson decay B Xs7. Very recently, 
^ \ just before the 60th birthday of Stefan Pokorski, this program has been formally completed. 
In the present paper, we summarize the existing results and discuss perspectives for further 
improvement of the accuracy of the Standard Model prediction for BR[i? -^sT]- 



Submitted to the special issue of Acta Physica Polonica B 
dedicated to Stefan Pokorski on the occasion of his 60th birthday. 



1 Introduction 



The radiative decay B X^'y is known to be extremely sensitive to the structure of fundamental 
interactions at the electroweak scale. It is dominantly generated by the Flavour Changing 
Neutral Current (FCNC) decay 6 — 57 that does not arise at the tree level in the Standard 
Model (SM). The leading order SM diagrams are shown in fig. |1|. 
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Figure 1: Leading-order Feynman diagrams for b 



57 in the SM. 



Many possible non-standard contributions (e.g., SUSY one- loop diagrams) are of the same 
order in electroweak interactions. They might remain important even for relatively heavy exotic 
particles. Consequently, 6 ^ 57 imposes severe constraints on extensions of the SM (see, for 
instance, ^-^). 

The inclusive branching ratio BR[5 — > Xs'~f] has been measured so far by CLEO 
BELLE and ALEPH 0]. The most accurate result is the one of CLEO, where photons 
with energies down to 2.0 GeV are included. Extrapolation towards lower photon energy cut- 
offs is performed following the phenomenological models of refs. [§, ^. 

When the photon energy cutoff is chosen to be 1.6 GeV in the 5-meson rest frame, the 
experimental world average read^ 



BR[B X,7 {E^ > 1.6 GeV)]exp = (3.12 ± 0.41) x 10" 



Within la, it matches the SM prediction [0, |T]_ 

BR[B Xs-f (E^ > 1.6 GeV)]sM = (3.57 ± 0.30) x 10"^ 



:i.2) 



One can see that the experimental and theoretical uncertainties are close in size. Without the 
inclusion of the Next-to-Leading Order (NLO) QCD corrections, the theoretical uncertainty 
in eq. (L2) would be around three times larger, and the constraints on new physics — much 
weaker. 

The program of the NLO calculation was outlined by Stefan Pokorski, Manfred Miinz and 
us in the article At that time, the only known results were the Leading Order (LO) ones 
that suffered from large scale uncertainties [10, |l|l. We analyzed these uncertainties in detail, 

^ By convention, contributions to B ~* Xsj from intermediate real ip and ?/;' are treated as background, 
while all the other cc states are included. 
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and enumerated calculations that still had to be done in the NLO case. Very recently, the last 
element of this NLO program has been completed fill]. In parallel to the QCD calculations, 
progress was being made in evaluation of the electroweak corrections, non-perturbative effects, 
as well as in collecting and analyzing the experimental data. 

In the present paper, we summarize all the contributions to the NLO QCD calculation 
of BR[5 —>■ Xg'j], and discuss perspectives for further improvement of the theoretical accu- 
racy. In particular, we point out the interplay between charm-quark mass uncertainties in the 
perturbative calculation and non-perturbative effects. 

Our article is organized as follows. The next section is devoted to a brief description of 
the history of perturbative calculations of QCD effects in 6 ^ s'j. In section |^, we summarize 
the electroweak corrections. Non-perturbative effects are discussed in section ^. The main 
theoretical uncertainties and possibilities for their elimination are the subject of section |^. 
Section contains our conclusions. 



2 The LO and NLO QCD calculations 

In a certain range of photon energy cutoffs, the width of the hadronic decay B —>■ Xgj is well 
approximated by the perturbative decay width 



(2.1) 



Arguments that support such a statement will be discussed in section 4. Until then, we shall 
restrict our discussion to the perturbative quantity ( |2.1|) . 

The framework for all the renormalization-group-improved perturbative analyses of 
b —>■ XP'^'"*°°7 is set by the effective Lagrangian 



eflf 



QCDxQED 



{u,d,s,c,b) + 



4G 



v2 i=i 



(2.2) 



It is obtained from the underlying theory (SM in our case) by decoupling of all the particles 
that are much heavier than the 6-quark. The Wilson coefficients Ci{^) play the role of coupling 
constants at the vertices Qi. The generic structure of the operators Qi is as follows: 

^SLa^-'bnF.,, I = 7, 



Q^ 



j^SLa^'^T-bnG 



a 



1,2, 



3,4,5,6, 



(g = u, d, s, c, b) 



(2.3) 



Here, Fj and F- denote various combinations of the colour and Dirac matrices (see, e.g., fl^ 



2 



The dots in eq. ( p.2| ) stand for UV-counterterms and non-physical operators that vanish 
by the QCDxQED equations of motion. In the present section, we neglect everything that 
is not important for 6 57 at the leading order in ttem, ^b/M]y, nis/mt, and Vut/Vct- This 
includes other operators Qi of dimension 5 and 6, higher- dimensional operators, as well as terms 
involving leptons. 

Let us assume that the decoupling of heavy particles is performed in the MS scheme, at the 
renormalization scale /zq ~ M^y. The values of Cj(/io) are found from the so-called matching 
conditions, i.e. by imposing equality of the effective- and underlying-theory Green functions 
at external momenta that are much smaller than masses of the decoupled particles. Next, 
the Wilson coefficients are evolved from fi = fiQ down to /i = yU;, ~ mf,, according to the 
Renormalization Group Equations (RGE) 

f^^Q{f^) = C,{fihM, (2.4) 

where the anomalous dimension matrix 7 is found from UV divergences in the effective theory. 
This procedure results in expressing the effective Lagrangian ( p.2| ) in terms of 

= Cr\,,) + )(^,) + (^)' + . . . , (2.5) 

where C\^\^h) depend on as only via the ratio r] = as(/io)/as(/i;,). Consequently, working at 
a fixed order in as, one truncates an expansion in powers of asifJ^b) rather than in powers of 
as{Mw) ln(M^/m^), as it would be the case without introduction of the effective theory. Thus, 
the behaviour of the perturbation series improves. This is the essence of the renormalizat ion- 
group improvement in the considered case. 



In the LO calculations, everything but Cf\^h) is neglected in eq. ( |2.5| ). At the NLO, 
one takes into account all the 0{as{^Jib)) contributions to T\b XP^''*°'^7], including those 
containing cl^\fib)- 

The Wilson coefficients encode information on the short-distance QCD effects due to hard 
gluon exchanges between the quark lines of the leading one- loop electroweak diagrams (fig. |l]). 
Such effects enhance the branching ratio BR[i? Xs'j] by roughly a factor of three, as first 
pointed out in refs. [T^ . 



A peculiar feature of the renormalization group analysis in 6 — 57 is that the mixing un- 
der infinite renormalization between the four-fermion operators Qi,...,Qg and the "magnetic 
penguin" operators Qt, Qs, which govern this decay, vanishes at the one-loop level. Conse- 
quently, in order to calculate the coefficients C-j{iib) and Cs{fJ^b) at LO, two-loop calculations 
are necessary. Such calculations were completed in ref. fl^. Earlier analyses [l^-|2y] contained 



additional approximations or were not fully correct. The results of ref. [|T^ were subsequently 
confirmed in refs. [^-[^]. 
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As pointed out in refs. |T^, |1]], the LO expression for r[b Xf^^^°^'y] suffers from large 
(~ ±25%) renormalization scale uncertainties. Therefore, matching the experimental accuracy 
of eq. (|1.1|) requires performing a complete NLO QCD calculation. This goal has been achieved 
in a joint effort of many groups: 

• Two-loop 0{as) corrections to the matching conditions C-j^hq) and C8(/io) were first 
calculated in ref. and subsequently confirmed by several groups p5[|-p^ 



Two-loop mixing and one- loop matching for the four-quark operators Qi,...,Qq were 
found in refs. [^-[^. In ref. I^S], these results were confirmed by recalculation in a 
different operator basis that is more suitable for 6 — > 57 analyses. 



Two-loop mixing in the sector (Qr^Qs) was calculated in ref. [3^|. These results have 
been recently confirmed p5 |. 



Three-loop mixing between the sectors (Qi, Qg) and (Qy, Qg) was evaluated in ref. p3| . 
It is currently being verified by another group 

The leading-order matrix elements {s'~fg\Qi\b) and the one-loop matrix element (s'jlQjlb) 
were calculated in refs. j^, 3^. Some of them were confirmed in ref. |37] where certain 



BLM corrections were included, too. 

Two-loop calculation of the matrix element {s'~f\Qi^2\b) was presented in ref. 



It has 



been recently verified and extended to the full basis of four-quark operators |]ri|, ^ . The 
one-loop matrix element (sjlQslb) has been found in refs. [^, |11|], too. 

It should be emphasized that all these ingredients enter not only the analysis oi B ^ Xg'y 
in the SM but are also necessary in extensions of this model. The corrections to the Wilson 
coefficients of the operators and Qs are also relevant for J3 — > XJ^l^ . 



3 Electroweak corrections 

The study of electroweak corrections begins with searching for terms that might be enhanced 
by large logarithms. Czarnecki and Marciano pointed out that large logarithms hi{ml/ml) 
are absent when is used in the overall normalization of F [6 ^ xp^''*°'^7] . 

Another type of large logarithm that might enhance some of the electroweak corrections is 
ln(m^/m^), i.e. the same logarithm that is responsible for the huge QCD enhancement of the 
6 — > S7 amplitude. Once [1 — as(/io)/as(/ife)] ~ 0.4 is treated as a quantity of order unity, the 
considered electroweak correction is formally of order (9(aem/cts), so it might be numerically 
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relevant, given the accuracy in eq. ( p..l| ). However, as demonstrated in refs. ^ through 

exphcit calculations, it turns out to be negligible (~ —0.7%). 

The articles |^ contain results for the complete electroweak corrections to the matching 
conditions Cj(yUo)- Some of them are proportional to aem(M^)/ sin^ 9^ ~ 0.034. Their effect 
on r[b Xf^^^°^'~f] amounts0to —1.5% for Mniggs = 115 GeV, and diminishes with increasing 
^Higgs- The authors of ref. |^ resolved the numerical discrepancy between refs. |^ and 
in favor of the latter. 

The only electroweak 0{aem) corrections that remain unknown at present are enhanced 
neither by large logarithms nor by 1/ sin^Oyy. Thus, we can be practically certain about their 
irrelevance. 



4 Non-perturbative effects 

The LO contribution to r[6 ^ xp^^"*°'^7] is given by the tree- level matrix element of the Qj op- 
erator. [| Let us temporarily assume that this operator is the only one in the effective Lagrangian 
(|2.2|) , and denote the corresponding contribution to the hadronic width hj T[B Xg'j]^'^'' oniy)_ 
In analogy to the analyses [0, ^ of the inclusive semileptonic decay B — > XuCu, one can 
apply the Operator Product Expansion (OPE) and Heavy Quark Effective Theory (HQET) to 
show that 



1 + ai — K + 02 — K + U \ 
mt mt \ mi 



(4.1) 



Here, Ai^2 ~ ^qcd standard HQET parameters. The value of A2 — 0.12 GeV^ is known 

from the measured B-B* mass difference. The value of Ai = -(0.27 ± 0.10 ± 0.04) GeV^ has 
been determined in ref. from the observed semileptonic -B-decay spectra (see ref. for 
more recent determinations). The coefficients oi and 02 can be calculated within perturbation 
theory,^ which yields |^, ^ 



1 9 
ai = - + 0{as{mi,)) and 02 = - - + 0{as{mb)). (4.2) 

The resulting 0{K^^^-Q/rn^) non-perturbative correction on the r.h.s. of eq. ( [4.1| ) amounts to 
around —3%. 

The relation ([4. 1|) still holds when a lower cutoff Eq is imposed on the photon energy in the 
i?-meson rest frame, provided Eq is not too close to the endpoint = '^2m_g^* — GeV. 
^ This number includes QED corrections to the matrix elements of Qi,2,7, too. 

^ In dimensional regularization, one-loop matrix elements of Q^.,...,Q(, may give LO contributions, too. 
However, they can be absorbed into the tree level matrix element of Qr with a redefined Wilson coefficient 
ailll- 
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The same refers to similar coefficients at higher orders in the (AQCD/m-6)-expansion. 



5 



1850801-007 




\ Data 

— Spectator Model 



1.5 2.5 3.5 4.5 

Ey (GeV) 

Figure 2: The B X^'-f photon spectrum observed by CLEO 



Acceptable values of Eq must correspond to much larger than Aqcd invariant masses of the 
recoiling hadronic state Xg. Fig. 3 in ref. suggests that Eq = 1.6 GeV is sufficiently low. 
More than 95% of the total r[6 — * Xf^'^^°'^'~f] originates from a peak that lays above such a 
cutoff. This peak is now clearly seen in the B X^'j spectrum observed by CLEO (fig. H). 
Its position corresponds to the photon energy in the leading two-body decay b sj. 

There is neither experimental nor theoretical need to consider photons below 1.6 GeV. 
They are practically unobservable at the inclusive level, because of the overwhelming b ^ c 
background. On the theoretical side, keeping not too small Eq facilitates the discussion of 
non-perturbative effects due to operators other than Qi- Of course, we have to admit that 
1.6 GeV is chosen arbitrarily. It could almost equivalently be 1.5 or 1.7 GeV. However, going 
up to the current CLEO cutoff of 2.0 GeV would increase uncertainties on the theoretical side. 
Data-driven extrapolation from the experimental cutoff to the theoretically preferred one is the 
right choice to make at present. 

The discussion of non-perturbative effects becomes much more complex when we take into 
account operators other than Qy. It is no longer possible to apply OPE in analogy to B ^ X^eu, 
because the 6-quark annihilation and the photon emission may now be separated in space-time 
by more than A~^^. 



The contribution of Qs to T[B X^^y] has been analyzed in ref. ||50| with the help of 
fragmentation functions. Important non-perturbative effects have been found for low E^ only, 
i.e. much below Eq = 1.6 GeV. Thus, with our cutoff, a reliable approximation is given by the 



Consequently, i?o-dependence in eq. (4.2) can be safely neglected. 
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perturbative contribution to r[6 — > Xf^^^°^'y] from the matrix elements of Qs- The accuracy of 
this approximation does not need to be known precisely, because the perturbative contribution 
of Qs is smaller than 3%. 

Similar conclusions can be drawn for the operators {srb){qr'q), where q = u,d,s. They 
are present inside Q3,...,Qq. Perturbative effects of their matrix elements are even smaller 
than that of Qs- As far as non-perturbative effects are concerned, one might worry about 
production of virtual vector mesons that convert to a real photon. However, creation of such 
transverse mesons is impossible in the factorization approximation because Q3,---,Qq contain 
no qcT^^q currents. Deviations from the factorization approximation are suppressed either by 



asirrii,) or by AQco/mf, |51]. This is sufficient to make them negligible here, given the smallness 
of |C'3,...,6(/^fe)| < 0.07, as compared to |Ci^2,7,8(Ai6)| — (0-5, 1, 0.3, 0.15). 

The operators Qs, ■■■,Qe contain {sTb)(bT'b) terms, too. The 6-quark loops are localized at 
distances much smaller A~^^ in space-time. Thus, they can undergo the same treatment as 
Qj, as far as non-perturbative effects are concerned. Since their perturbative contributions are 
minor, the non-perturbative ones are totally negligible. 

Charm quark loops are the most difficult to analyze. Factorization is not sufficient here 
because 2mc/mh is not a small number. Moreover, non-factorizable contributions may be 
numerically important because the Wilson coefficients Ci and C2 are not small at all. 

Let us begin with tracing down possible contributions from intermediate real cc states. Our 
cutoff £^0 = 1-6 GeV implies that the invariant mass of the final Xg state is smaller than 
rriri^ + rriK- Consequently, real cc states might occur only before the photon emission, i.e. in a 
cascade decay: B YccX^^^ followed by Ycc — > X^'^^j. 

The importance of such processes can be tested in the case Fee = i^, because separate 
experimental data on both (inclusive) components of the cascade decay are available. For low 
Eq, the resulting branching ratio of the intermediate ip contribution is larger than the one 
in eq. (p^). It gets reduced to (a few)xlO~^ for Eq = 1.6 GeV, and becomes negligible for 
Eq = 2.0 GeV m . 




Figure 3: Charm loop contribution to 6 — > s'jg. 



The models used by CLEO p, |5^ to extrapolate from 2.0 GeV to lower cutoffs do not 
include the intermediate ip contribution. They are based on perturbative calculations, in which 
the only diagram (fig. ^ that might correspond to this contribution affects r[6 — > XP'^'"*°°7] 
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by 1.7% only (for Eq = 1.6 GeV). Consequently, the procedure applied by CLEO is consistent 
with treating the intermediate ip contribution as background. 

Identical arguments work for Higher cc states might produce higher energy photons. 
However, radiative charm annihilation processes in all the cc states except and ip' have negli- 
gible branching ratios. Thus, it does not really matter whether we consider their contributions 
as background or not. Whatever decision is made, its effect is expected to be less than the 
1.7% perturbative contribution from the diagram in fig. ^ 

Having discussed the real intermediate cc states, we proceed to the virtual ones. Nei- 
ther infrared nor collinear singularities occur in the perturbative contributions of cc loops to 
r[b —>■ Xf^^^°^'-f] at NLO. Thus, according to the common wisdom, one expects that these per- 
turbative results give reasonable estimates to the corresponding contributions to r[i? — > X^'j], 
up to corrections of order O {Aqq^ / rricfi) ■ 




Figure 4: One-loop matrix element that vanishes for the on-shell photon. 

The actual situation is somewhat more complicated, because the leading one-loop diagram 
(fig. ^) vanishes for the on-shell photon. However, it becomes non-vanishing when a soft gluon 
is attached to the c-quark loop. Such a gluon may originate from the decaying B meson. 
Thus, one finds a non-perturbative effect ^ that is not approximated in any sense by the 
corresponding perturbative null. Fortunately, it can be expressed within HQET in terms of a 
series 



m 



in which the n > 1 terms are likely to be negligible, because the coefficients 6„ decrease rapidly 
with n 0. The calculable leading 0{A'^^^/ml) term enhances the decay width by around 



2.5% |5l 

The perturbative 0{as) results described in section ^ include non-vanishing two-loop dia- 
grams with cc loops, e.g. the ones obtained by adding a virtual gluon to the diagram in fig. |^. 
The corresponding non-perturbative effects are expected to be suppressed by both as{mb) and 
Aqch/ ^c,b- Thus, at the first glance, they might seem irrelevant. However, it remains an open 
question whether their suppression is numerically sufficient. No quantitative estimates of such 
non-perturbative effects have been performed so far. We shall discuss this issue in more detail 
at the end of the next section. 
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5 Phenomenological discussion 



In the present section, we shall discuss the two main uncertainties in the present-day SM 
prediction for B — > X^j. The analysis of ref. |^ will be largely followed. 
The prediction (|1.2|) is obtained from the formula 



= BR[B X^eu 
in which the following substitutions are made 




V[B 



Xucv] 



th 



(5.1) 



^ V[B ^ X^eu] 



T[6-^XP^''*°°7 {E^>Eq)\ 



th 



' Y\B ^ Xuev\ \ 

y\b^xm)^^ 



Jj^^partoiigj^l 



NLO 



non-pert urbative 
corrections (4.3) 



NNLO 



known (9(A2/'n2^) 
corrections 



,(5.2) 



(5.3) 



Such ratios are introduced in order to minimize uncertainties in eq. ( ^.1|) that originate from 
the CKM angles and the overall factors of m^. The use of 6 — > m transitions is motivated by the 
fact that eq. (|5.3|) is known at the NNLO, while convergence of the perturbation series and non- 
perturbative effects are more easily controlled in eq. ( |5.2| ) than in r[i? ^ Xs^\lY\B — >■ XcCz/]. 
The 0(AQpj3/m^) terms from eq. ( |4.1D have canceled in the ratio (|5^ ) with the analogous 
corrections to r[5 — > XuCp]. On the other hand, once the quark masses are expressed in terms 



of the hadronic ones, the ratio (|5.3| ) depends on both Ai and A2. 




h s b s b s b s 

Figure 5: Charm loop contributions to the matrix elements of four-quark operators. 

The main uncertainty in the perturbative ratio on the r.h.s. of eq. ( ^.2|) originates from the 
two-loop diagrams with charm quarks presented in fig. |^. Such diagrams are the only source of 
mc-dependence of this ratio. Since the higher-order (NNLO) QCD corrections are unknown, the 
renormalization scheme for ttIc remains arbitrary, at least within a certain class of "reasonable" 
schemes that do not artificially enhance the unknown corrections. As argued in ref. 0, the 
uncertainty in eq. (p..2|) stemming from this scheme- dependence can be accounted for by setting 
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mc/rrib = mc^jJif^^ / rn\^ in the two-loop diagrams,^ and varying the scale /x between rric and 
nil,. Such a variation is the dominant source of the error in eq. ( p. .21) . 

One could remove the considered uncertainty by calculating three-loop diagrams obtainable 
from fig. ^ by adding one more virtual gluon. UV-divergent parts of such diagrams have been 
already found in the process of calculating the NLO anomalous dimensions ||2^. Evaluating 
the finite parts would constitute an extremely tedious task, though not totally impossible, if 
numerical integration was applied. Finding the remaining NNLO corrections would be relatively 
simpler, given that fully automatized analytical methods are now available [^, pl[] . 

However, before undertaking such an ambitious task, one should make sure that all the 
non-perturbative effects are really under control. The main worry are the doubly-suppressed 
corrections mentioned in the last paragraph of section ^. So far, they have been neither esti- 
mated nor included in the theoretical error. They are related to precisely the same two-loop 
diagrams with charm quarks (fig. |^). Numerical importance of non-local parts of those dia- 
gram^ can be illustrated by the fact that the r.h.s. of eq. (|1 . 2|) changes by 35% when rric is 
shifted from the original value of 0.22 m;, to the threshold for charm pair production nic = ^mf,. 
A AqcD/^^fe-suppressed non-perturbative effect on the top of such a large perturbative contri- 
bution might not be negligible. Unfortunately, no systematic methods have yet been developed 
for calculating corrections of this type. 



6 Conclusions 

In the present paper, we have summarized the existing calculations of perturbative and non- 
perturbative contributions to the inclusive weak radiative B meson decay. We have pointed out 
that both the main perturbative uncertainty and the most worrisome non-perturbative effects 
have their origin in the fact that non-local charm quark loop contributions are particularly large. 
Removing the perturbative uncertainty due to mc-dependence would be extremely tedious, but 
not totally impossible. However, developing a method for systematically estimating the related 
non-perturbative effects is desirable in advance. 

The present agreement at the ~ 10% level between the experimental ( |1.1| ) and theoretical 
(|1.2|) determinations of BR[i? Xs'~f] implies that clear signatures of new physics in this 
observable are not likely to be found in the foreseeable future. The importance of improving the 
accuracy on both the experimental and theoretical sides follows from the need for strengthening 
the b —>■ s'-f constraints on beyond-SM theories. Such constraints are likely to be crucial in 
identifying the origin of new physics effects that we expect to encounter in the LHC era. 

^ Here, ml^ stands for the 6-quark mass in the so-called "IS-scheme" |Q. It is defined as half of the 
perturbative contribution to the T mass. 

^ By non-local we mean those parts that cannot be removed off-shell by finite local counterterms. 



10 



7 Acknowledgements 



We would like to thank K. Baranowski, C. Bobetli, K.G. Chetyrkin, A. Czarnecki, P. Gambino, 
M. Jamin, A. Khodjamirian, A. Kwiatkowski, M.E. Lautenbacher, M. Miinz, S. Pokorski, 
N. Pott, J. Urban and P.H. Weisz for fruitful collaboration in the calculations that were relevant 
for determining the SM prediction for b —>■ sj. A. J.B. is grateful to the Max Planck Institute for 
Physics for financial support during the Planck02 Symposium in Kazimierz. His work presented 
here has been supported in part by the German Bundesministerium fiir Bildung und Forschung 
under the contract 05HT1WOA3 and the DFG Project Bu. 706/1-1. M.M. was supported in 
part by the Polish Committee for Scientific Research under the grant 2 P03B 121 20. 

References 



A.J. Buras, M. Misiak, M. Miinz and S. Pokorski, Nucl. Phys. B 424 (1994) 374 fEep 
ph/931 134511 . 



M. Misiak, S. Pokorski and J. Rosiek, In the review volume "Heavy Flavours II", eds, 
A.J. Buras and M. Lindner, Adv. Ser. Direct. High Energy Phys. 15 (1998) 795 fhep- 
ph/9703442| 



P. Gambino and M. Misiak, Nucl. Phys. B 611 (2001) 338 ||hep-ph/0104034 



J.R. Elhs, D.V. Nanopoulos and K.A. Olive, Phys. Lett. B 508 (2001) 65 ||hep-ph/0l025M | 



S. Chen et al. (CLEO Collaboration), Phys. Rev. Lett. 87 (2001) 251807 [|hep-ex/0108032 |. 



K. Abe et al. (Belle Collaboration), Phys. Lett. B 511 (2001) 151 ||hep-ex/0 103042 



R. Barate et al. (ALEPH Collaboration), Phys. Lett. B 429 (1998) 169. 

A. Ah and C. Gre ub, Z. Phys. C 49 (1991) 431, Phys. Lett. B 259 (1991) 182, Phys. Lett. 
B 361 (1995) 146 ||hep-ph/950637l . 



A.L. Kagan and M. Neubert, Eur. Phys. J. C 7 (1999) 5 [ |hep-ph/98053Q3 
A. Ah and C. Greub, Z. Phys. C 60 (1993) 433. 



A.J. Buras, A. Czarnecki, M. Misiak and J. Urban, Nucl. Phys. B 631 (2002) 219 fEep 
ph/0203135| 



S. Bertolini, F. Borzumati and A. Masiero, Phys. Rev. Lett. 59 (1987) 180. 
N.G. Deshpande et a/., Phys. Rev. Lett. 59 (1987) 183. 



M. Ciuchi ni et al, Phys. L ett. B 316 (1993) 127 [ |hep-ph/9307SBl , Nucl. Phys. B 421 
(1994) 41 ||hep-ph/9311357|] . 



11 



B. Grinstein, R.P. Springer and M. B. Wise, Phys. Lett. B 202 (1988) 138, Nucl. Phys. B 
339 (1990) 269. 

R. Grigjanis, P.J. O'Donnell, M. Sutherland and H. Navelet, Phys. Lett. B 213 (1988) 
355, Phys. Lett. B 286 (1992) 413 (E). 

G. Cella, G. Curci, G. Ricciardi and A. Vicere, Phys. Lett. B 248 (1990) 181. 
M. Misiak, Phys. Lett. B 269 (1991) 161. 



K. Adel and Y.P. Yao, Mod. Phys. Lett. A 8 (1993) 1679 f|hep-ph/9302244 

M. Misiak, NucL Phys. B 393 (1993) 23. 

M. Misiak, NucL Phys. B 439 (1995) 461 (E). 



G. Cella, G . Curci, G. Ricciardi and A. Vicere, Phys. Let t. B 325 (1994) 227 jE^ 
ph/9401254i , Nucl. Phys. B 431 (1994) 417 ||hep-ph/9406"20^] . 



K.G. Chetyrkin, M. Misiak and M. M iinz, Phys. Lett. B 400 (1997) 206, Phys. Lett. B 
425 (1998) 414 (E) ||hep-ph/96123Tl , Nucl. Phys. B 518 (1998) 473 [|liep-ph/9711266 . 



K. Adel and Y.P. Yao, Phys. Rev. D 49 (1994) 4945 ||hep-ph/9308"349 |. 

C. Greub and T. Hurth, Phys. Rev. D 56 (1997) 2934 ||hep-ph/9703"349|l . 

A.J. Buras, A. Kwiatkowski and N. Pott, Nucl. Phys. B 517 (1998) 353 ||liep-ph/9710336[ . 



M. Ciuchini, G. Degrassi, P. Gambino and G.F. Giudice, Nucl. Phys. B 527 (1998) 21 
[|hep-ph/9710335 . 



C. Bobeth, M. Misiak and J. Urban, Nucl. Ph ys. B 567 (2000) 153 [|hep-ph/9904413|] , 
Nucl. Phys. B 574 (2000) 291 [|hep-ph/991022q] . 



G. Aharelli, G. Curci, G. Martinelh and S. Petrarca, Nucl. Phys. B 187 (1981) 461. 
A.J. Buras and P.H. Weisz, Nucl. Phys. B 333 (1990) 66. 

A.J. Buras, M. Jamin, M.E. Lautenbacher and P.H. Weisz, Nucl. Ph ys. B 370 (1992) 69 
(Addendum - ibtd. B 375 (1992) 501), Nucl. Phys. B 400 (1993) 37 |hep-ph/ 92 11301 . 



M. Ciuchini, E. Franco, G. Martinelh and L . Reina, Phys. L ett. B 301 (1993) 263 |Eep 
ph/9212203i , Nucl. Phys. B 415 (1994) 403 ||hep-ph/9304"25^] . 



K.G. Chetyrkin, M. Misiak and M. Miinz, Nucl. Phys. B 520 (1998) 279 ||liep-ph/9711"280| . 
M. Misiak and M. Miinz, Phys. Lett. B 344 (1995) 308 ||hep-ph/940lM3^ . 
P. Gambino, M. Gorbahn and U. Haisch, in preparation. 
N. Pott, Phys. Rev. D 54 (1996) 938 ||hep-ph/9512"252| . 



12 



Z. Ligeti, M.E. L uke, A.V. Manohar and M.B. Wise, Phys. Rev. D 60 (1999) 034019 
||hep-ph/ 9903305 1 . 



C. Greub, T. Hurth an d D. Wyler, Phys. Lett. B 380 (1996) 385 ||hep-ph/9602281|] , Phys. 
Rev. D 54 (1996) 3350 ||hep-ph/9603"40l . 



A. J. Buras, A. Czarnecki, M. Misiak and J. Urban, NucL Phys. B 611 (2001) 488 |Eep| 
ph/0105160| 



A. Czarnecki and W.J. Marciano, Phys. Rev. Lett. 81 (1998) 277 ||hep-ph/9804252 



K. Baranowski and M. Misiak, Phys. Lett. B 483 (2000) 410 [|hep-ph/9907427 



P. G ambino and U. H aisch, JHEP 0009 (2000) 001 ||hep-ph/0007259|l , JHEP 0110 (2001) 
020 ||hep-ph/010905g] . 



A. Strumia, Nucl. Phys. B 532 (1998) 28 ||hep-ph/9804"274l] . 

J. Chay, H. Georgi and B. Grinstein, Phys. Lett. B 247 (1990) 399. 

A.V. Manohar and M.B. Wise, Phys. Rev. D 49 (1994) 1310 ||hep-ph/9308"24^ . 

A.H. Hoang, Z. Ligeti and A.V. Manohar, Phys. Rev. Lett. 82 (1999) 277 Eep 



ph/9809423|| , Phys. Rev. D 59 (1999) 074017 ||hep-ph/ 98 11239 



C. W. Bauer and M. Trott, |hep-ph/0205"039| . 
A.F. Falk, M.E. Luke and M.J. Savage, Phys. Rev. D 49 (1994) 3367 ||hep-ph/930g^ 



LL Bigi et al, DPF92 Proceedings (Batavia, November 1992) ||hep-ph/9212227[ . 

A. Kapustin, Z. Ligeti and H.D. Pohtzer, Phys. Lett. B 357 (1995) 653 [|hep-ph/9507248|l . 



M. B eneke, G. Buchal la, M. Neubert and C. T. Sach rajda, Phys. Rev . Lett. 83 (1999) 
1914 ||hep-ph/9905"3ll] , NucL Phys. B 591 (2000) 313 ||hep-ph/ 0006121 . 



A. Khodjamirian and M. Misiak, unpubhshed (2000). 
E. H. Thorndike (CLEO collaboration), |hep-ex/0206067 . 
M.B. Voloshin, Phys. Lett. B 397 (1997) 275 ||hep-ph/9612"483 



A. Khodjam irian, R. Riickl, G. StoU and D. Wyler, Phys. Lett. B 402 (1997) 167 jEi^ 
ph/9702318| 



Z. Ligeti, L. Randall and M.B. Wise, Phys. Lett. B 402 (1997) 178 [ |hep-ph/9702"322l . 

A.K. Grant, A.G . Morgan, S. Nussinov and R.D. Peccei, Phys. Rev. D 56 (1997) 3151 
||hep-ph/ 97023801 . 



[58] G. Buchalla, G. Isidori and S.J. Rey, Nucl. Phys. B 511 (1998) 594 ||hep-ph/ 970525^ 



13 



[59] T. van Ritbergen, J. A. Vermaseren and S. A. Larin, Phys. Lett. B 400 (1997) 379 fhep- 
I ph/970139D| , Phys. Lett. B 405 (1997) 327 ||hep-ph/970328l . 



[60] K.G. Chetyrkin, Acta Phys. Polon. B 28 (1997) 725 ||hep-ph/ 96 10531 



[61] M. Steinhauser, Comput. Phys. Commun. 134 (2001) 335 ||hep-ph/0009ll^ |. 



14 



